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A general overview of the measurement of long-range multiplicity correlations measured 
by the STAR experiment in Au+Au collisions at RHIC is presented. The presence of 
^ long-range correlations can provide insight into the early stages, and the type of matter 

, produced in, these collisions. These measurements have been made in Au+Au collisions 

at v / sjvjv = 200 and 62.4 GeV. These results indicate a relatively large long-range corre- 
' lation is produced in Au+Au collisions compared to a pp baseline at y/sjvjv = 200 GeV. 

r ^ i A weaker long-range correlation is seen as a function of incident energy. Further, corn- 

ed , parison of the onset of the long-range correlation to the calculated percolation density 

parameter at y/s mn = 200 GeV is presented. 
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s_i ■ 1. Introduction 



The study of particle correlations has been suggested as a mea ns to search for a 
possible phase transition in ultra-relativistic heavy ion collisions JH2E2 These parti- 
cle correlations have the potential to probe the early stages of heavy ion collisions. 
If a deconfined phase of quarks and gluons exists in these collisions, the presence 
of these partonic degrees of freedom could directly influence the results of a corre- 
lation measurement. One such measurement is the study of long-range multiplicity 
correlations that extend over a range in pseudorapidity (77). A linear relationship 
has been experimentally shown in hadron-hadron collisions relating the multiplicity 
in the forward direction (Nf) to the average multiplicity in the backward direction 

<N b (N f )>=a + bN f (1) 



The coefficient b is referred to as the correlation coefficient. The value, a, is 
a measure of uncorrelated particles. In the definition provided in "Eq. ([1])" the 
maximum value of b is one, corresponding to complete correlation between the 
produced particles. For 6 = 0, all particles are uncorrelated. This relationship can 
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be expressed in terms of the following expectation valued 

< N f N b >-<N f ><N b > Djf 

<N*>-<N f > 2 D) f {> 

where D 2 j and D 2 j are the backward-forward and forward-forward dispersions, re- 
spectively. This result is exact and model independent.^ All quantities in "Eq. ([2])" 
are directly measurable. The large pseudorapidity and angular coverage of the 
STAR experiment at RHIC provides an ideal means to measure particle correlations 
with precision.^ 

In phenomenological particle production models, the hadronization of color 
strings formed by partons in the collision account for the produced soft particles. 
In the Dual Parton Model (DPM) , at low energies these strings are formed between 
valence quarksP At higher energies, there are additional strings formed due to the 
increasing contribution from sea quarks and gluons. In the DPM, this production 
mechanism is predicted to produce long-range multiplicity correlations in pseu- 
dorapidity. A long-range correlation could be an indication of multiple elementary 
inelastic collisions in the initial stages of an A+A collision.^ Additionally, the Color 
Glass Condensate (CGC) model predicts long-range correlations.® The CGC model 
is a description of high density gluonic matter at small-x. In the CGC model, the 
incoming, Lorentz contracted nuclei are envisioned as two sheets of colliding colored 
glassP In the initial stage of the collision, the color fields are oriented transverse 
to the beam direction. After the collision these color fields realign longitudinally. 
This produced state is referred to as the Glasma, an intermediate state between 
the CGC and Quark Gluon Plasma (QGP)P^ These longitudinal fields provide the 
origin of the long-range correlation and are similar to the strings in the DPM. 



2. Data Analysis 

The data utilized for this analysis are from year 2001 (Run II) ^smn = 200 GeV 
Au+Au and year 2004 ^Jsnn = 62.4 GeV collisions at the Relativistic Heavy 
Ion Collider (RHIC), as measured by the STAR (Solenoidal Tracker at RHIC) 
experiment.^ The main tracking detector at STAR is the Time Projection Cham- 
ber (TPC) The TPC is located inside a solenoidal magnet generating a constant, 
longitudinal magnetic field. For this analysis, data was acquired at the maximum 
field strength of 0.5 T. All charged particles in the TPC pseudorapidity range 
0.0 < \rj\ < 1.0 and with px > 0.15 GeV/c were considered. This 77 range was 
subdivided into forward and backward measurement intervals of total width 0.2rj. 
The forward-backward correlations were measured symmetrically about rj = with 
varying 77 gaps (measured from the center of each bin), A?y = 0.2, 0.4, 0.6, 0.8, 1.0, 
1.2, 1.4, 1.6, and 1.8 units in pseudorapidity. Fig Q] depicts a schematic representa- 
tion of the forward-backward correlation and its measurement. The collision events 
were part of the minimum bias dataset. The minimum bias collision centrality was 
determined by an offline cut on the TPC charged particle multiplicity within the 
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Rapidity interval 

Fig. 1. Schematic diagram of the measurement of a forward- backward correlation. 

range \r/\ < 0.5. The centralities used in this analysis account for 0-10, 10-20, 20-30, 
30-40, 40-50, 50-60, 60-70, and 70-80% of the total hadronic cross section. An addi- 
tional offline cut on the longitudinal position of the collision vertex (v z ) restricted 
it to within ±30 cm from z = (center of the TPC). Approximately two million 
minimum bias events at y/sjfjj 1 = 200 GeV and 0.7 million events at ^/snn = 
62.4 GeV satisfied these requirements and were used for this analysis. Corrections 
for detector geometric acceptance and tracking efficiency were carried out using a 
Monte Carlo event generator and propagating the simulated particles through a 
GEANT representation of the STAR detector geometry. 

In order to eliminate the effect of statistical impact parameter (centrality) fluctu- 
ations on this measurement, each relevant quantity ((Nf), (Nb), (Nf) , and (NfNb)) 
was obtained on an event-by-event basis as a function of STAR reference multiplic- 
ity, N c h- A linear fit to (Nf) and (Nb), or a second order polynomial fit to (Nf) and 
(NfNb) was used to extract these quantities as functions of N c h- Tracking efficiency 
and acceptance corrections were applied to each event. These were then used to cal- 
culate the backward-forward and forward-forward dispersions, D\f and D"jf, binned 
according to the STAR centrality definitions and normalized by the total number 
of events in each bin. An autocorrelation exists if the centrality definition over- 
laps with the measurement interval in r/ space. Therefore, there are three r\ regions 
used to determine reference multiplicity: < 0.5 (for measurements of A77 < 1.0), 
0.5 < |?/| < 1.0 (for measurements of A77 > 1.0), or | -^7 1 < 0.3 + 0.6 < |?/| < 0.8 (for 
measurement of A77 = 1.0). This is possible due to the flat r\ yield of the STAR 
TPC within the measurement ranges being considered. 
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3. Results 

Fig. [2] presents the mean multiplicity (Nf) and (Nb), in a forward and backward 
pseudorapidity region along with (Nf) and (NfNb) as lunctions of reference mul- 
tiplicity, N c h- The results are calculated event-by-event for every unit N c h. The 
(Nf) and (Nb) demonstrate a linear dependence over much of the centrality range. 
Only statistical uncertainties arc present and are less than 1%, except for the high- 
est multiplicities. The results in Fig. [2] are uncorrected for tracking efficiency and 
acceptance. Due to statistical limitations, it is not possible to apply corrections for 
every value of N c h- One value for the correction, calculated for each centrality, is 
applied to the resultant values (Nf) , (Nb) , (Nf) , and (NfNb) for every event in 
that centrality. Therefore, all events falling within a particular centrality have the 
same correction. 

Fig.[3]shows the results for the correlation coefficient b and the backward-forward 
and forward-forward dispersions, D\* and D 2 ^, for most central (0-10%) Au+Au 
collisions at ^/snn = 200 GeV. The correlation is approximately flat as a function 
of the pseudorapidity gap. It is expected that if only short-range correlations are 
present (mostly from sources such as cluster formation, jets, or resonance decay), 
then b would decrease rapidly as a function of A77EI It has been argued that the 
presence of these long-range correlations at larger than Ai] = 1.0 cannot be due to 
final state effects, and are therefore due to the conditions in the initial statePl This 



a Note that there are methods to estimate the short-range contribution and remove it, thereby 
providing an absolute measure of the "pure" long-range component. This is beyond the scope of 
the current presentation, but is discussed in the following!^ 



February 8, 2008 21:57 WSPC/INSTRUCTION FILE ws-ijmpe 



Energy Dependence of Short and Long-Range Multiplicity Correlations in Au+Au Collisions from STAR 5 



0.8 
0.7 
0.6 
0.5 
1 0.4 
0.3 
0.2 
0.1 


450 
400 



£ 350 

o 

g» 

(]) 250 

Q. 

</) 200 



150 
100 
50 




(a) 



(b) 



0-10% Au+Au 



STAR Preliminary 



u bf 



STAR Preliminary 



1 1.5 

A T) 



Fig. 3. (a) Correlation coefficient, b and (b) Dispersions, and as a function of the 

pseudorapidity gap (A77) in 0-10% most central y'sjvlv = 200 GeV Au+Au collisions. If only 
short-range correlations were present, the expectation is a rapidly falling b value with increasing 
An. The approximately flat behavior of b indicates the presence of a long-range correlation. 



may be phenomcnologically similar to fluctuations larger than the event horizon 
being created during an inflationary phase in the early universe that later become 
smaller than the event horizon!^ 

A comparison of the correlation coefficient, b, in 0-10% most central ^snn — 
200 and 62.4 GeV is shown in Fig. [U There is a slight difference in the evolution of 
the correlation at small values of At], but for A77 > 0.6 the value of the correlation 
coefficient is approximately flat as a function of Arj for the two energies. However, 
for the ytsNN = 62.4 GeV data, the value of b is lower by up to w 15%. It is expected 
that b has a dependence on incident energy (and atomic number) of the colliding 
nuclei.^ Data for Cu+Cu collisions at the same energies are currently under study. 

The measurement of a long-range correlation demonstrates a coherent phe- 



February 8, 2008 21:57 WSPC/INSTRUCTION FILE ws-ijmpe 



6 Terence Tarnowsky 



x: 0.8 
c 

j= 0.6 

</) 

■S 0.4 
re 

0) 

5 0.2 
O 



-0.2 



STAR Preliminary 



Au+Au 200 GeV 



Au+Au 62.4 GeV 



0.2 0.4 0.6 0.8 



1 1.2 

A T| 



1.4 1.6 1.8 



Fig. 4. (a) Comparison of the correlation coefficient, b, as a function of the pseudorapidity gap 
(Ajj) in 0-10% most central ^/sjvjv = 200 and 62.4 GeV Au+Au collisions. The lower energy 
data exhibits a similar plateau as the 200 GeV result, but at a lower absolute number. This could 
indicate a smaller long-range correlation in Wsmn = 62.4 GeV Au+Au collisions. 



nomenon in the initial state of these heavy ion collisions. It may then be possible to 
correlate this effect with that of another collective phenomenon, percolation. The 
particle production mechanism in high-energy heavy ion collisions can be des crib ed 
by the hadronization of color strings stretched between the colliding nucleifSl At 
high energies it is possible for partons from different nucleons to begin to overlap 
in the p lane tr ansverse to the propagation direction of the nuclei, forming clusters 
(Fig. [5|)! 15 | 16| ^ the percolation threshold, a cluster that spans the entire system 
is formed. This point can be characterized b y the criti cal percolation density pa- 
rameter, r] c . r\ c has been predicted theoreticall y^ 7 | 18 J19| an ^ percolation density 
parameter itself h as bee n measured experimentally in heavy ion collisions by the 
STAR experiment EOEH jt has been shown by these measurements that the criti- 
cal value, r] c , for the "local" prediction {r] c = 1.72) is exceeded in most centralities 
of Au+Au collisions at ^snn = 200, while for y / 7/\w = 62.4 Au+Au collisions 
all centralities lie at or below the critical percolation threshold. A slightly lower 
value for rj c is predicted in the infinite (continuum) system. In that case, central 
Au+Au at ytsNN = 62.4 GeV will lie above r\ c . It is thus important to try and 
experimentally constrain r\ c . If the presence of long-range correlations and perco- 
lation are related, it may be possible to constrain r\ c by examining the long-range 
correlation as a function of the measured percolation density parameter, shown in 
Fig. El The b shown here is the long-range component onlyU^I it is seen that the 
extracted long-range component of the correlation coefficient b crosses zero in the 
area around a percolation density parameter value of two. This is higher than most 
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Fig. 6. Extracted long-range correlation, &Li}Ci as a function of the percolation density parameter, 
r], for Au+Au at y/sjvjv = 200 GeV. At the point where the long-range component fades, rj is ~ 
2. Most theoretical predictions for eta c lie below this. 



theoretically predicted values for r\ c . Though this result is intriguing, further study 
is required. 

4. Summary 

Results have been presented on the energy dependence of the forward-backward 
correlation coefficient, b, in Au+Au collisions. It is seen that the value of b is ap- 
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proximately flat over a wide range in pseudorapidity gap, Arj, at ^/snn = 200 
GeV. This indicates a long-range correlation as a function of Arj. Were only short- 
range correlations present, a rapid decay in the value of b would be expected as 
A77 increases. These long-range correlations can be ascribed to multiple elementary 
inelastic collisions. These are predicted in models such as the Dual Parton Model 
and the Color Glass Condensate/Glasma phenomenology. The long-range correla- 
tion appears to persist at t/snn — 62.4 GeV, though at a lower value than ^/snn 
= 200 GeV. Further investigations are proceeding to determine an experimental 
range for the critical percolation density parameter, r\ c . At the point where the 
extracted long-range correlation value goes to zero, the value of the percolation 
density parameter is w 2. This is larger than several theoretical predictions for r/ c . 
Additional study is required to characterize this relationship. 
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